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I. Program of Research for the Second Quarter 

During this past quarter: (1) the heliocentric velocity dependencies of 

the OH photodissociation lifetime and the (0-0) band fluorescence rate have 
been incorporated into the cometary radical model, (2) preliminary evaluation 
of the effects of radiation pressure and the heliocentric distance dependent 
parent velocity of CN production rates deduced from photometry has been 
performed, and (3) study of radical scale length data has been continued. 

1. OH Photochemistry 

Schleicher and A'Hearn (1982, 1984) have calculated the dependencies of 
both the OH photodissociation lifetime and the OH solar fluorescence rate on 
heliocentric velocity. At a distance of 1 AU from the sun, the OH lifetime 
can vary from <J2 x 10^ s to -4 x 10^ s depending on the heliocentric velo- 
city. The (0-0) band fluorescence rate varies by more than a factor of 6 for 
different heliocentric velocities. Figures 1 and 2 show these dependencies 
which have been incorporated in to the cometary radical particle-trajectory 
model. 


2. CN Production Rate 

The cometary radical model, which was updated last quarter to include the 

- 1/2 

heliocentric distance dependent parent molecule velocity, v = 0.58 r^ 

(Delserame 1982), has been applied to recalculate the production rates cf CN in 

Comet West (A'Hearn et al. 1977, A'Hearn and Cowan 1980). A'Hearn et al. 

(1977) had originally used an r^ + ^ law for a Haser model parent scale length 

for CN and found that the production rate for CN varied with heliocentric 

distance as Combi and Delsemme (1980b) then measured CN parent scale 

lengths in Comets Bennett and West and found a dependence of r H * * which 

o 

they concluded was reasonably consistent with an r^ law. Both they and 

2 

A'Hearn and Cowan (1980) then adopted the r^ law and found a CN production 
rate varying approcimately as r H ~ n where n * 1.6±0.3. 

For this study, the CN production rates have been calculated using three 
different model descriptions and compared with the Haser model parameters 
A'Hearn and Cowan have adopted. The results of these four cases are shown in 
Table l. The first case is a Haser model, using the Average Random Walk Model 
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EXCITATION RATE FOR OH (0-0) 
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HELIOCENTRIC VELOCITY (KM/S) 



Table 1 




CN Production in Comet 

West 





Model 

Q! (s' 

l } (a) 

(b) 

n 


A 

Hearn and Cowan Haser Model 

1.1 X 

10 27 

-1.61 

1) 


Average Random Walk Model 

1.4 x 

10 27 

-2.00 

2) 

Monte 

Carlo Particle Trajectory Model 
(no radiation pressure) 

1.2 x 

10 27 

-1.92 

3) 

Monte 

Carlo Particle Trajectory Mode’ 
(with radiation pressure) 

1.3 x 

10 27 

-1.87 


Production Rate at I AU 
Exponent in power law Q = Qj r 
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(Combi and De’ _mme 1980a) corrections to the scale lengths. Here the r H 
variation in parent velocity and the isotropic ejection of the daughter 
radicals is taken into account in an approximate way. The second case is a 
Monte Carlo Particle Trajectory Model (MCPTM) which explicitly calculates the 
effects included in the first case. And the third case is a MCPTM which also 
includes radiation pressure (Combi 1980). In all cases the heliocentric 
velocity dependence of the CN(O-O) band excitation (Tatum and Gillespie 1977) 
is included. 

The largest effect is that of the variation of the parent velocity which 
steepens the slope of the resulting production law and is present in all three 
cases. The explicit inclusion of isotropic ejection (case 2) and radiation 
pressure (case 3) each flatten the slope somewhat but are nonetheless impor- 
tant. A comparison of the original Haser model calculations and the full 
MCPTM (Case 3) for the CN production in Comet West is shown in Figure 3. 

3. C 2 and CN Scale Lengths 

Although Haser model scale lengths cannot be used to give direct physical 
quantities (lifetimes), they can be useful tools to characterize the apparent 
size of the source and decay regions of cometary radicals. Newburn and 
Spinrad ( 1984) have presented Haser scale lengths for the parents of C 2 , C-j 
and CN in several periodic comets. These were obtained by two-point filter 
photometry: one point centered on the photometric nucleus and another 

displaced some known distance from the nucleus (typically 1 to 4 x 10^ km). 
Since Haser's model is basically a three parameter model (the production rate 
and two scale lengths), Newburn and Spinrad had to assume the values of 
A'Hearn (1982) for the radical decay length. From two points they could then 
calculate both a parent scale length and the production rate. Although errors 
could arise from both this assumption and that of spherical symmetry (i.e., no 
radiation pressure distortion), a large number of observations may average out 
some of the random discrepancies. 

Figure 4 shows a plot of the CN parent scale length versus heliocentric 
distance for the combined data sets of Newburn and Spinrad (1984), Combi and 
Delsemme (1980) and Delsemme and Combi (1983). There is no gross systematic 
difference between the two data sets. The scale lengths of Newburn and 
Spinrad were determined from the two-point photometry of smaller short period 
comets at medium to large heliocentric distances, whereas our data were 
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Haser Model Scale Length Ckm) 


C 


T 1 1 1 1 — TT 



Helioce 


a P/ENCKE 
A P/TUTTLE 
O P/STEPHAN-OTERMA 
▲ BRADFIELL) 

■ P/SWIFT-GEHRELS 
x BENNETT 

+ KOHOUTEK 
• WEST 





determined from the average of the sunward and antlsunward brightness profiles 
of larger new or very long period comets at small to medium heliocentric 
distances. Taking all the data together, a power law can be fit which yields 
the straight line in Figure 4 and which has the form X p * 1.6 x 10^ r^*^ for 
X p in km and r^ in AU. 

The slope for this fit has an uncertainty of about ±0.3 generally 

12 15 

excluding both the typically adopted r H or r H laws. An r H * aw » on the 

n 

other hand, would be expected for a photodissociation lifetime (x * r H ) a 
variable parent velocity (v « r H ”®*^). 

Figure 5 shows the variation of the C 2 parent scale length with helio- 
centric distance for the combined data sets. Unlike the case for CN, though, 
there does seem to be a systematic difference. The scale lengths of Newburn 
and Spinrad tend to be both larger and exhibit a flatter slope than do those 
of Combi and Delsemme (1985). Taken together, a power law in r^ can be fitted 
to all the data which have the form X * 2.1 x 10^ r^**"*. ® ur data taken 

A 1ft 

alone has the form X^ = 1.4 x 10 r u i * , whereas the data of Newburn and 

p „ H ’ 

Spinrad yield X p - 3.8 x 10 A r H °*'. 

At this point, there could be two reasons for such a difference. First, 
the C 2 spatial distribution in the smaller short period comets may be substan- 
tially different from that in the more active new and very long period comets. 
If C 2 were produced primarily by gas phase chemical reactions or an icy grain 
source as suggested by A' Hearn and Cowan (1980), there should certainly be 
differences between the two populations of comets. Second, there could be 
differences between the two model fitting procedures which did not surface in 
the CN data. The C 2 radical experiences nearly twice the acceleration and is 
likely moving at a lower speed than is CN (Combi and Delsemme 1985). Since 
the two-point photometry method does not average out the asymmetry, systematic 
differences could result. 
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II. Prograa of Research for the Next (^tarter 


Research activ'.cies In the next quarter will concentrate on (1) develop- 
ment of a multiple collision algorithm for the MCPTM to treat neutral-neutral 
collisions in the inner coma, and (2) preliminary model runs for the OH 
distribution. 
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Figure Captions 


Figure 1. Photodissociation Lifetime of OH as a function of Heliocentric 
Velocity (Schleicher and A'Hearn 1982). 

Figure 2. Excitation Rate for the OH(O-O) Band. The heliocentric velocity 
dependence is shown for heliocentric distances of 0.25 AU and 1.00 AU 
(Schleicher and A'Hearn 1984). 

Figure 3. CN Production Rate vs. Heliocentric Distance. The CN Production 
rates for Comet West are shown as computed hoth with the Haser model 
scale lengths of A'Hearn and Cowan (1980) [0] and with the Monte Carlo 
particle trajectory model which includer radiation pressure and a vari- 
able parent velocity [+]. Intermediate cases are discussed in the text. 

Figure 4. CN Parent Scale Length vs. Heliocentric Distance. The points 

plotted for Comets P/Encke, P/Tuttle, P/Stephan-Oterma and Bradfield were 
computed from two-point filter photometry by Newburn and Spinrad 
(1984). Those for Comets Bennett, Kohoutek and West were computed from 
entire brightness profiles (Combi and Delsemme 1980b, Delsemme and Combi 
1983). The solid line is the best fit power law to the combined data 
sets. 

Figure 5. C 2 Parent Scale Length vs. Heliocentric distance. The points 
plotted for Comets F/Encke, P/Tuttle, P/Ste?han-Oterma and P/Swift- 
Gehrels were computed from two-point filter photometry by Newburn and 
Spinrad (1984). Those for Comets Burnham, Tago-Sato-Kosaka, Bennett and 
Kohoutek were computed from entire brightness profiles (Combi and 
Delsemme 1985). The solid line is the best fit power law to the combined 
data sets. 
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